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ABSTRACT 

High-resolution atmospheric flow simulations of the tidally-locked extrasolar giant planet, HD 209458b, 
show large-scale spatio-temporal variability. This is in contrast to the simple, permanent day/night (i.e., 
hot/cold) picture. The planet's global circulation is characterized by a polar vortex in motion around 
each pole and a banded structure corresponding to ^3 broad zonal (east-west) jets. For very strong 
jets, the circulation-induced temperature difference between moving hot and cold regions can reach up 
to ^1000 K, suggesting that atmospheric variability could be observed in the planet's spectral and 
photometric signatures. 

subject headings: planetary systems - planets and satellites: general - stars: atmospheres - turbulence 



1. INTRODUCTION 

Roughly 100 gaseous giant planets are currently known 
to orbit nearby sun-like stars 9 . Many of those planets 
are located at very small orbital distances from their par- 
ent stars, where tidal forces are thought to maintain a 
rotation rate synchronous with the orbit — thus producing 
permanent day and night sides on the planet. This sit- 
uation presents a new regime of atmospheric circulation, 
not encountered in our solar system: slowly-rotating giant 
planets, which are continuously exposed to intense stellar 
heating on the same side. Measuring the resulting temper- 
ature structure on these planets is a major goal of current 
and future observational programs. In this Letter, we re- 
port on high-resolution, fully-turbulent global simulations 
of the atmospheric flow on HD 209458 b, presently the only 
close-in extrasolar giant planet (CEGP) with a measured 
mass (M p ) and radius (R p ). 

The parent star of HD 209458 b shows discernible bright- 
ness decrements every 3.5 days, due to occultations by the 
planet as it transits across the star. This property has 
recently led to precise measurements of M p and R p (Char- 
bonneau et al. 2000; Henry et al. 2000; Mazeh et al. 
2000), confirming its giant nature. It has also allowed the 
detection of sodium absorption, providing the first probe 
of the planet's atmosphere (Charbonneau et al. 2002). 
According to the standard planetary formation picture, 
HD 209458 b is expected to have formed at a large dis- 
tance (> 1 AU) from its parent star (see, e.g., Boss 1996) 
and migrated inward (Goldreich & Tremaine 1980; Lin et 
al. 1996; Murray et al. 1998), quickly (<; 10 Myr) reach- 
ing its present distance of only 0.046 AU from the star (see 



e.g., Burrows et al. 2000). There, it was forced by tidal 
effects to permanently present the same face to its star 
(see, e.g., Goldreich & Soter 1966) — as the Moon does to 
the Earth. From this synchronization, the rotation period 
of the CEGP is known (same as its orbital period of 3.5 
days). However, unlike our Moon with its insignificant 
atmosphere, HD 209458 b is expected to possess vigorous 
meteorology and associated horizontal transport of heat 
and chemical species, due to the presence of a thin, stable 
(radiative) atmospheric region above the convective inte- 
rior (Guillot et al. 1996; Seager & Sasselov 1998). 

2. MODEL 

We model the stable region of HD 209458 b's atmosphere 
as a shallow layer of hydrostatically balanced, frictionless 
gas — moving under the influence of gravitational and Cori- 
olis accelerations. The motion of such a layer enveloping 
a planet is governed by the equivalent barotropic formu- 
lation of the shallow-water equations on a rotating sphere 
(Salby 1989): 
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where v = v (A, ip, t) is the horizontal velocity with A and 
if the longitude and latitude, respectively; h = h(X, ip, i) 
is the thickness of the modeled layer, proportional to the 
temperature; f — 2U,sanp is the Coriolis parameter, where 
£7 is the rotation rate of the planet; k is the local unit vec- 
tor normal to the planetary surface; g is the gravitational 
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acceleration; T a = —gVi](l — exp{— t/r a }) cos(<p) cos(A) 
and Td = ~{h — hE)lT ( i are, respectively, the adiabatic and 
diabatic thermal forcing, representing uneven hemispheric 
heating and cooling due to synchronization; r a and are 
characteristic e-folding times; Ke = H + rj cos(ip) cos(A) is 
the equilibrium thickness; 77, and H are constants; and, 
K, = R/c p , where R is the specific gas constant and c p is 
the specific heat at constant pressure. 

The nonlinear equations (1) and (2) describe circula- 
tions in which the motion field is vertically aligned over 
one or more pressure scale heights, H p . They are column- 
integrated representation of the equations used in gen- 
eral circulation models and are solved numerically us- 
ing the highly-accurate pseudospcctral algorithm (Eliasen, 
Mechenhauer & Rasmussen 1970; Orszag 1970). Approx- 
imately 150 simulation runs with up to T341 (1024x512 
grid) resolution have been performed to explore the full 
physical and numerical parameter-space available for the 
HD 209458 b atmosphere. The details of the full explo- 
ration are described elsewhere (Cho et al. 2003). Here, 
we emphasize only the main robust features from the 
large number of simulations — namely, the emergence of 
long timescale variability due to moving polar vortices and 
"thermal" spots of high contrast from the background tem- 
perature field. 

High-resolution, shallow-layer barotropic models in 
spherical geometry have been successfully used to model 
the atmospheric dynamics of the solar system gaseous gi- 
ant planets (Williams 1978; Cho & Polvani 1996a), as well 
as stratospheric (stable region) phenomena on the Earth 
(see, e.g., Juckes & Mclntyre 1987). An alternate model 
of observed cloud level circulation on Jupiter as a surface 
expression of deep convective columns, oriented parallel to 
the planetary rotation axis, also exists (Sun et al. 1993; 
Schubert & Zhang 2000). However, only the shallow-layer 
model has thus far qualitatively reproduced the essential 
features (bands, zonal winds, persistent spots, and anticy- 
clonic dominance) on all of the giant planets in our solar 
system using only the observed values of physical param- 
eters (Cho & Polvani 1996a). 

For HD 209458 b, we adopt parameter values where 
known (Charbonneau et al. 2000; Henry et al. 2000; 
Mazeh et al. 2000): /C = 0.29, p8m s" 2 , R p = 10 8 m, 
n = 2.1 x 10" 5 s" 1 , and H = 2 x 10 6 m. The last 
(= H p /JC) represents the global mean layer thickness, 
whose value corresponds to the equivalent blackbody tem- 
perature (Tatm ~ 1400 K for a Bond albedo, ^4 = 0.2; the 
results presented here, however, are insensitive to the pre- 
cise choice of r atm and A). The values for two required 
parameters — 77, the amplitude of hemispheric thermal forc- 
ing, and U, the characteristic wind speed related to the 
global mean kinetic energy — are not known and hence 
varied. From run to run, rj is varied from to 0AH, 
corresponding respectively to no heating and maximum 
substellar-point heating in the absence of atmospheric mo- 
tion; U is varied from 50 m s _1 to 1000 m s , roughly the 
observed value for Jupiter (Ingersoll 1990) and the value 
expected from a simple thermal wind balance (Showman 
& Guillot 2002), respectively. Note that the latter value 
leads to a nonlinearly-balanced wind field which can lo- 
cally approach the sound speed (^2700 m s _1 ); on all ob- 
served atmospheres in our solar system, U does not exceed 
~400 m s" 1 (Bcatty & Chaiken 1990). We have carefully 



studied the effects of U on the variability and thermal 
contrast and find that, while a variability in the flow is 
a generic feature for HD 209458 b, a significant thermal 
contrast generally requires a small 77 or a large U. In the 
following, we give an example of the expected generic flow 
pattern and then describe a typical flow-induced thermal 
variability which may be observationally relevant. 

3. RESULTS 

Fig. 1 shows an instantaneous atmospheric flow field on 
HD 209458 b at dynamic equilibrium from one of our sim- 
ulations with U = 400 m s _1 and 77 corresponding to a 
maximum day/night temperature difference of 180 K. The 
flow exhibits considerable complexity and variability: it is 
characterized by a broad equatorial band containing thin, 
undulating filamentary structures at low and middle lat- 
itudes and displaced, translating circumpolar vortices at 
the rotation poles. The large undulations at low latitudes, 
which can roll up into small vortices, are due to high- 
amplitude Rossby (planetary) waves. In addition to being 
responsible for weather, the planetary waves efficiently stir 
and mix tracers (e.g., clouds, chemical species, potential 
vorticity) when they break, homogenizing the region in 
which the breakings occur (Juckes & Mclntyre 1987). In 
contrast, polar vortices sequester air inside, their sharp 
boundaries serving as robust barriers to material transport 
between their interiors and the well-mixed interstitial re- 
gions. The overall flow is similar to that of Earth's winter 
stratosphere. 

The complex flow of Fig. 1 corresponds to a global cir- 
culation pattern of 3 zonal jets, shown in Fig. 2. Jets are 
atmospheric flow structures, which are obtained by lon- 
gitudinally averaging the eastward wind at each latitude. 
The pattern in Fig. 2 naturally evolves from an initially 
random flow without jets, with the basic profile forming 
approximately at day 20. The emergence of jets from a 
random field is a generic feature of turbulent flows in a dif- 
ferentially rotating shallow layer (Rhines 1975; Williams 
1978; Cho & Polvani 1996b). For the extreme values of 
U considered with all other parameters held fixed, 2 or 
4 jets are also possible. However, once formed, the jets 
do not change qualitatively — even under applied thermal 
forcing — for up to several hundred HD 209458 b days, the 
duration of the simulation. 

The presence of a low number of broad jets, along 
with polar vortices, leads to a circulation pattern which 
is markedly different from that of cloud-top Jupiter. Al- 
though CEGPs are sometimes dubbed "hot Jupiters" , the 
circulation pattern is more like that of Uranus or Nep- 
tune, or even Earth, than that of Jupiter or Saturn, which 
have ~10 narrow jets (Ingersoll 1990). Dynamically, this 
is expected given roughly the factors of ten greater T atm 
and ten slower rotation rate of HD 209458 b compared to 
Jupiter. The combined effect, giving the flows in Figs. 1 
and 2, is characterized by the Rossby deformation scale, 
Ld= (ICgH) 1 / 2 1 /, which is to be compared with R p . Gen- 
erally, Ld/ 'R p <C 1 is required for a Jupiter- like multiple- 
jet circulation pattern without pronounced polar vortices 
at the cloud-top. For HD 209458 b, L D /Rp~l. 

Jets and vortices/eddies can strongly influence the over- 
all thermal structure by transporting heat. Since the state 
of the HD 209458 b's cloud-top is not likely to be at rest 
and driven solely by stellar heating, we delay the applica- 



CHO ET AL. 



3 



tion of thermal forcing (in contrast to Showman & Guil- 
lot 2002) until the flow field self-organizes into the pat- 
tern in Fig. 2. Dynamically, this is a likely flow state 
of HD 209458 b in the presence of weak or uniform heat- 
ing (i.e., at pre-synchronization). The forcing is slowly 
increased to the specified value of r\ over a characteris- 
tic radiative equilibrium time, r r = 10 days, appropriate 
for ~1 bar pressure level on HD 209458 b. At this depth, 
T r ^> t c (where t c = R p /U is the characteristic circulation 
time), signifying a large distance advection of airmass be- 
fore it warms or cools appreciably. The results presented 
here do not change qualitatively with the value of r r (~ 
0-100 days), using either an adiabatic {T a ,T a = r r ) or a 
diabatic (J r d,Td = T r ) representation of thermal forcing. 

At each point on the planet, the layer thickness is related 
to the temperature field, T(A, <p, t) = hg/c p , via the hypso- 
metric relation (Holton 1992; Salby 1989). In our simula- 
tions with small r\ or large U, the T field exhibits a generic 
feature which is in direct contrast to the simple perma- 
nent day/night picture. Fig. 3 shows the T distribution 
(from the same run of Figs. 1 and 2) at several different 
times. Note that the warmest and the coldest regions are 
not located at the sub-stellar and the anti-stellar "thermal 
poles" (at the equator), respectively. Instead, the temper- 
ature extrema are located near the rotation poles, inside 
a coherent hot/cold spot pair at each pole. Nor are the 
extrema stationary: they revolve around the poles with a 
period of ^25 days in this case. The movement leads to a 
distribution that for a time has the temperature minimum 
actually on the "hot" day side and the maximum on the 
"cold" night side. In addition, the temperature difference 
between the two extrema is large and strongly asymmet- 
ric; here, the difference is ^300 K, with the diffuse hot 
region ~70 K above and the sharp cold region ^230 K be- 
low T atm . For the case with U — 1000 m s _1 , the difference 
reaches ^1000 K, with a temperature minimum of ^800 K 
below T atm at the core of the cold spots. 

4. DISCUSSION 

The presence of high-contrast hot and cold spots on 
HD 209458 b induces spatio-temporal variability, which 
may produce detectable fluctuations in observational sig- 
natures. Sensitive enough infrared flux measurements 
(e.g., such as those possible with SIRTF) could reveal vari- 
ability in time during an orbit of the planet as different 
faces are seen from Earth, as well as from orbit to orbit as 



the spots (polar vortices) revolve about the rotation poles. 
The lack of detectable variability would thus point toward 
cither an obscuring uniform haze overlying the modeled 
region or an inefficient conversion of stellar irradiation to 
atmospheric kinetic energy. In the latter case, the formed 
spots are weak (i.e., small temperature/thickness pertur- 
bation) and thermal forcing which produces bulging in the 
modeled layer of more than several percent overwhelms 
any temperature variability due to atmospheric motion. 
The equilibrium day/night temperature is then robustly 
maintained. 

The extreme conditions inside the spots in our calcula- 
tions suggest several additional potential observables since 
absorption levels, albedo, intrinsic thermal emission, and 
presence, type, or height of clouds could all be different 
within the different spots. The spatially integrated spec- 
trum, therefore, can be different from the uniform planet 
case. For example, enhanced levels of CH4/CO abundance 
inside the cold spots are possible, given the value of T atm 
used (Seager et al. 2000). Similarly, condensates, such 
as MgSiC>3 (enstatite), may also be found inside the dif- 
fuse hot spots, where temperatures may be high enough 
(Sudarsky et al. 2000; Seager et al. 2000). If the temper- 
ature at radiative equilibrium is actually higher than as- 
sumed (e.g., ^2000 K), the opposite situation may occur — 
MgSiC>3 may be found in the cold spots. In addition, by 
sequestering chemically active species and periodically ex- 
posing them to the stellar irradiation (as the spots revolve 
around the poles) , the spots could also affect atomic num- 
ber densities and could be part of the explanation for the 
recently observed low abundance of Na I on HD 209458 b 
(Charbonneau et al. 2002). 
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Fig. 1. — Two views of the dynamical flow tracer, potential vorticity (Holton 1992), at day (=year) 55 from our T341 (1024x512 grid 
resolution) simulation of the atmospheric circulation of HD 209458 b: a) orthographic projection centered at the anti-stellar point (AS) on the 
night side and b) polar-stereographic projection centered at the north pole (NP). 1 PVU = 4xl0~ 27 s -1 !!! 1 '^'. The global flow is characterized 
by two circumpolar cyclonic (rotating in the same direction as the planet — counter-clockwise in the figure) vortices at high latitudes and 
high-amplitude planetary waves at low latitudes. 
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Fig. 2. — Steady circulation profile from the flow of Fig. 1: longitudinally-averaged zonal (eastward) wind velocity. Zonal jets emerge from 
a random initial condition due to conversion of small-scale stirring/eddies to planetary waves in the atmosphere (as seen in Fig. 1), a generic 
feature of shallow-layer turbulence under differential (latitudinally-dependent) rotation /. This basic profile has emerged approximately at 
day 20 of the simulation and does not change for the duration of the simulation (several hundred HD 209458 b days). On HD 209458 b, the 
jets are fewer, broader, and stronger than those observed on Jupiter. This is because HD 209458 b is hotter and rotating more slowly than 
Jupiter. Hemispheric thermal forcing is applied on this robust flow. 



Fig. 3. — Four successive night-side views of the temperature field in orthographic projection (on HD 209458 b days indicated). A persistent 
"thermal dipole", associated with the circumpolar vortices, revolves around the pole and provides high thermal contrasts that may be 
detectable in the planet's observational signatures. In contrast to the simple, permanent day/night picture, the hottest atmospheric region 
is, at times, on the night side of the planet. These spots are long-lived, asymmetric (cold regions are more well-defined and have larger 
amplitude), and can serve as large areas of atmosphere with distinct chemistry and thermodynamics. The global equilibrium blackbody 
temperature of the atmosphere is taken to be 1400 K. 
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